Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN KTpX style file vl.4) 



Spectral evolution of magnetic flares and time lags in accreting 
black hole sources 



Juri Poutanen ' and Andrew C. Fabian 

1 Stockholm Observatory, S-133 36 Saltsjobaden, Sweden 

2 Institute of Astronomy, Madingley Road, Cambridge, CB3 OHA 



Accepted 1999 April 28. In the original form 1998 November 23 



ABSTRACT 

We present a model for the short timescale spectral variability of accreting black holes. It 
describes well the time-averaged spectra as well as temporal characteristics such as power- 
density spectrum, time/phase lags, and coherence function of Cygnus X-l. We assume that 
X/7-rays are produced in compact magnetic flares at radii ^ 100GM/c 2 from the central 
black hole. The tendency for magnetic loops to inflate and detach from the underlying ac- 
cretion disc causes the spectrum of a flare to evolve from soft to hard due to the decrease 
of the feedback from the cold disc, so causing time delays between hard and soft photons. 
We identify the observed time lags with the evolution timescales of the flares, which are of 
the order of the Keplerian timescale. We model the overall temporal variability using a pulse 
avalanche model in which each flare has a certain probability to trigger a neighbouring flare, 
thus occasionally producing long avalanches. The duration of the avalanches determines the 
Fourier frequencies where most of the power emerges. 

Key words: accretion, accretion discs - radiation mechanisms: non-thermal - stars: flare - 
stars: individual (Cygnus X-l) - X-ray: general - X-ray: stars 



1 INTRODUCTION 

The X-ray emission of many sources harbouring a black hole 
shows strong variability, on time-scales of milliseconds to hours 
for objects like Cyg X-l and other Galactic black holes (GBHs), 
and minutes to days for an active galaxy like MCG-6-30-15. The 
time-averaged spectra are often interpreted as the result of the 
Comptonization of soft photons in a hot static electron cloud (e.g. 
Shapiro, Lightman & Eardley 1976; Sunyaev & Triimper 1979). 
Surprisingly, even modern X/7-ray data of one of the best studied 
bright sources, Cyg X-l, are relatively well explained in terms of 
such models (Gierlinski et al. 1997; Zdziarski et al. 1998; Poutanen 
1998), despite the observed rapid variability. Detailed studies of the 
X-ray colours, however, have shown rapid spectral changes (Nolan 
et al. 1981; Negoro et al. 1994), which imply the presence of rapid 
changes in the physical conditions of the source. 

The properties of the observed variability are somewhat puz- 
zling. In the case of Cygnus X-l, the power-density spectrum 
(PDS) shows that most of the power emerges at Fourier frequen- 
cies corresponding to time-scales between 0.1 - 5 s (Belloni & 
Hasinger 1990; van der Klis 1995). Studies of the autocorrelation 
function revealed also the existence of strong correlations for lags 
up to ~ 5 s (Nolan et al. 1981; Meekins et al. 1984). These time- 
scales exceed by orders of magnitude the characteristic time scales 
of the processes operating close to the black hole. 
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Another intriguing feature is the Fourier-frequency-dependent 
time lags between hard and soft photons (Miyamoto et al. 1988; 
Miyamoto et al. 1991; Cui et al. 1997; Nowak et al. 1999). Since it 
is generally believed that Comptonization is responsible for spec- 
tral formation, these lags are sometimes attributed to the delays ex- 
pected from scattering in a hot cloud. Hard photons are the result 
of more scattering and so emerge after, or lag behind, softer ones 
(Payne 1980). In order to explain lags of order ~ 0.1s, a huge 
hot cloud of size R ~ 10° - 10 10 cm is then required. Comp- 
tonization in a uniform cloud, however, produces lags independent 
of Fourier frequency (Miyamoto et al. 1988). This led Kazanas, 
Hua & Titarchuk (1997) to propose an inhomogeneous density dis- 
tribution of the hot cloud. While solving one problem, however, 
this model raises another, more serious, one of how physically to 
support a cloud of such huge size for which most of the energy is 
at the outer boundary while all the dissipation is at the centre. Be- 
sides, when modelling the time lags, it is also often assumed that 
the Compton cloud is static and it is the soft photon input which 
varies, contrary to what is observed (Miyamoto et al. 1991). 

Miyamoto & Kitamoto (1989) were first to propose the spec- 
tral evolution as a solution for these problems. However, in their 
model an ad hoc scaling between time-scales of hard and soft pulses 
is assumed. Here we present a new model for the spectral variabil- 
ity of accreting black holes and assume that the basic unit of energy 
release is a flare which defines the shortest variability time-scale. 
The flares are assumed to be produced from an active region in a 
stochastic manner which we model as a pulse avalanche. Individual 
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flares are assumed to move away from the disc, so that the fraction 
of the reprocessed soft flux that crosses the flare emission region 
(= ER) decreases with time during the life of a flare. This produces 
the spectral variability that causes the time delays between hard and 
soft photons of the order of the flare time-scale. 



2 PHYSICAL PICTURE 

Magnetic flares have long been considered as the source of the 
rapid variability of Cyg X- 1 and other similar accreting black holes 
(Galeev, Rosner & Vaiana 1979; Haardt, Maraschi & Ghisellini 
1994; Stern et al. 1995; Di Matteo 1998; Beloborodov 1999). Mag- 
netic fields are amplified by the differential motion in the disc, rise 
up into the corona where they reconnect and liberate their energy in 
a flare. Exactly how this process takes place is not clear, but there 
are indications that they expand and detach from the disc (see Mi- 
kic, Barnes & Schnack 1988; Amari, Luciani & Tagger 1996; Ro- 
manova et al. 1998). One expects the overall evolution time-scale 
of the magnetic field configuration to be of the order of Keplerian 
time-scale (Romanova et al. 1998). Since most of the energy in the 
accretion disc is dissipated within radii JS 50R g from the central 
black hole (here R g = 2GM/c 2 ), these time-scales lie between 
~ 1 ms and a few tenths of a second for a M ~ IOMq black hole. 

The characteristic size of the magnetic loop (giving an upper 
limit for the size of the ER) is of the order of a turbulent cell size, 
equal to the half-thickness of the accretion disc (Galeev et al. 1979). 
For the innermost part of the accretion disc it is then R ~ rhR g 
(Shakura & Sunyaev 1973), where rh = Mc 2 / Xsdd ~ 1 in the 
case of Cyg X-l, and can be smaller if a large fraction of energy is 
dissipated above the disc, when the disc height is reduced (Svens- 
son & Zdziarski 1994). 

The reconnection of the magnetic field lines in one loca- 
tion can force the neighbouring field configurations into an unsta- 
ble state causing subsequent reconnections and, possibly, a long 
avalanche. The evolution of an individual avalanche may roughly 
parallel that of a Solar active region (Lu & Hamilton 1991). 



3 SPECTRAL EVOLUTION OF A MAGNETIC FLARE 

We follow here the scenario where magnetic fields are moving 
away from the underlying accretion disc in the course of the flare 
and estimate how the flare spectrum varies in time. We assume that 
a heating rate lh(H) (I = Lot / (-Rm e c 3 ) is the compactness pa- 
rameter, and L is the luminosity) is a function of height, H (mea- 
sured in units of the characteristic size of the ER, R), and is zero 
both when the ER is at the surface of the disc and when it is far 
away, reaching maximum in between. 

The hot electrons produced by the dissipation are cooled by 
Comptonization of intrinsic soft photons produced in the disc by 
viscous dissipation (with corresponding compactness h,o) and pho- 
tons produced by reprocessing in the disc of flare radiation. In the 
beginning of the flare, when lh is smaller than Z Sj o, the soft photons 
are just the average background of luminosity l St o, and the spec- 
trum emitted by the flare is very soft. When l h is large, the hard 
X-rays produced in the ER heat up the cold disc in the vicinity of 
the flare, so creating extra soft photons which then dominate the 
cooling. The emitted spectrum defined by the geometry of the sys- 
tem is hard (Stern et al. 1995; Svensson 1996). The luminosity of 
the soft photons that crosses the ER is l s (H) = l s ,o + D(H)lh(H). 



The height (and time) dependent feedback factor, D(H), deter- 
mines what fraction of the hard luminosity after reprocessing in 
the disc returns to the ER. During the whole feedback phase when 
h(H) S> h,o the emission region moves away from disc thus de- 
creasing D(H). This causes a continuous increase of the Comp- 
ton amplification factor, A = l h (H)/l s (H), and a corresponding 
hardening of the spectrum. In the very last phase, when intrinsic 
disc photons dominate cooling, the spectrum becomes soft again. 

A hardening of the spectrum during the phase of maximal en- 
ergy dissipation agrees well with the observation of short time-scale 
spectral evolution of Cyg X-l (Negoro et al. 1994). The luminosity 
in hard photons reaches a maximum later than the soft luminosity 
does. In order to show that the model described above can qualita- 
tively be applied to GBHs, we have computed the temporal evolu- 
tion of the flare spectrum and the time-averaged spectrum using the 
code of Coppi (1992) and compared it with that of Cyg X-l. For 
an ER of R R g , we can assume that the spectrum is produced 
instantaneously (i.e. there are no time delays due to Compton up- 
scattering). This is a good approximation as long as characteristic 
time-scale of the flare r 3> R/c. The spectral evolution then can be 
described as a sequence of steady states. At any given moment the 
spectrum is defined by three parameters: l s , lh, and Thomson opti- 
cal depth of the ER, tt . The Compton reflection component is also 
assumed to be produced in the vicinity of the flare, thus varying 
together with the Comptonized radiation. 

The feedback factor can be parametrized as D(H) = 
Do/(l + 3H 2 ), where the uncertainties in geometry of the ER, 
structure of the cold disc and the angular distribution of the flare 
radiation are incorporated into the factor Do- The temperature of 
the disc that produces soft seed photons for Comptonization also 
depends on the distance of the ER from the disc: 

CTT-KcrsD7r y l + SH^J 

and equals to T S (H) = 0.25[Z S , + r)lh(H)/(l + 3H 2 )] 1/4 keV 
for R ~ 3 • 10 6 cm (~ R g , for a IOMq black hole). The factor 
rj combines the uncertainties in the angular distribution of the hard 
radiation and albedo for the reflection from the disc. 

Assuming l s ,o = 0.01, the temperature of the cold disc be- 
tween flares is T Sj o = 80 eV. We fix the compactness of a flare 
at the peak of the energy dissipation at /a, max = 10 (Haardt et al. 
1994; Beloborodov 1999). As long as the ratio lh,ma.^/h,o is large, 
the spectrum depends weakly on our choice of Z^.max (and Z s ,o), but 
depends strongly on the feedback and the height where most of the 
energy dissipation occurs. We further assume, for illustration, that 
the heating rate varies with time as lh{t) oc (t/r)^ exp(— tjr), and 
the ER moves away from the disc with a constant velocity, so that 
H — Hot/r. The exact shape of the flare profile (value of £) does 
not affect any of the qualitative results. Hard spectra observed in 
GBHs require most dissipation to occur at a height approximately 
equal to the size of the ER (Svensson 1996) if the X-ray emis- 
sion is isotropic (i.e. rj ~ 1/2). Thus we take £ = 2 (ad hoc), 
H = 1/| = 1/2, t] = 1/2, and fit D . Altogether, we have 
only two free parameters, Do, and the Thomson optical depth, tt, 
that determine the flare spectrum. To be specific, Do determines the 
slope of the spectrum, while tt determines the electron tempera- 
ture and the spectral cutoff. We note that parameters Ho and Do 
are strongly correlated, i.e. if one decreases Ho, the same spectrum 
can be produced decreasing Do- In the case of the thermal plasma 
(assumed here) and I ~ 10, the contribution from the e ± pairs is 
negligible for the spectrum observed in Cyg X-l. 
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Figure 1. (a) Spectral evolution of a magnetic flare. Time resolved spectra (without Compton reflection) are presented by dotted curves; marks are times 
in units of r from the beginning of the flare. Time-averaged Comptonized spectrum plotted by a thin solid curve corresponds to the Compton amplification 
A ~ 10. The Compton reflection is computed by convolving the flare spectrum with the angle-dependent Green functions of Magdziarz & Zdziarski (1995) 
with inclination fixed at 30°. Time-averaged spectrum of Cygnus X-l (simultaneous Ginga and OSSE data from June 1991, data set 1 in Gierlihski et al. 1997) 
is plotted by crosses and the best fit with the flare model by a solid curve (x 2 /dof = 50.0/75). Interstellar absorption is removed when plotting the model 
spectrum. The best fit model parameters, Do = 0.50+q'q?, and tt = 1.8^q'^. Here for spectral fits we use XSPEC vlO (Arnaud 1996), and the errors are 
given for a 90 per cent confidence interval (A% 2 = 2.7). The amount of Compton reflection Q/(2tt) = 0.301q'q3, the EWof Fe line is 100 eV, and absorber 
column density N H = (0.9+° ?,) ■ 10 22 cm~ 2 . The contribution to the \ 2 is shown at the bottom panel, (b) The flare light curves (in arbitrary units) at 3, 
9, 27, and 81 keV are presented by dotted curves; the heating rate, J/,(t), by solid curve. Temperature of the ER, kT e (dashed curve), peaks at the time when 
Compton amplification A(t) is largest. Pulses at different energies are almost self-similar, i.e. can be represented by the same (t/r^;) 2 exp(— */tb) law with 
time-constants at photon energy E approximately following the logarithmic dependence, te ~ t[1 + i \n(E / E v )], where r is the time-constant of the 
energy dissipation £/,(t) and E p is the photon energy where EFg peaks. Thus the delay between photons of energies E and Eq is At j \n(E/Eo). 



With these assumptions, we compute the spectral evolution of 
the flare (see Fig. la) and find that the time-averaged spectrum of 
the flare gives a very good fit to the data of Cyg X-l. If there were 
anisotropy breaks in the time-resolved spectra (Stern et al. 1995) 
due to the anisotropy of the soft photon injection into the ER (ne- 
glected here), the spectral evolution would smooth them out. In the 
isotropic case, the time-averaged Comptonized spectrum is slightly 
concave and generally cannot be described by a Comptonized spec- 
trum with a fixed optical depth and temperature. 

From the spectral evolution we can compute the flare light 
curves in various energy bands (Fig. lb). Due to the soft-to-hard 
evolution, soft photons come first, while the hard photon flux (at ~ 
100 keV) follows the evolution of the heating rate and peaks later. 
This delay of hard photons relative to soft ones causes time lags to 
appear in the composite light curves from an ensemble of flares. 



4 THE PULSE AVALANCHE MODEL 

GBHs as well as many other complex physical systems often show 
1// noise in their PDS (Press 1978). There is no accepted theory 
for the origin of this noise. One possibility is that a system develops 
into a self-organized critical state (Bak, Tang & Wiesenfeld 1987; 
Mineshige, Ouchi & Nishimori 1994; Mineshige & Negoro 1999). 
In general, the correlations between signals (flares) of very differ- 
ent time scales are required. We model such correlations using a 
version of a simple stochastic pulse avalanche model due to Stern 
& Svensson (1996) originally invented to describe temporal char- 
acteristics of gamma ray bursts. 

We assume that some flares occur spontaneously according to 
a Poisson distribution with a mean rate of A flares per second (as 



in the shot noise model, Terrel 1972). All flares can stimulate the 
production of flares, the number following a Poisson distribution 
of mean number p. Each spontaneous flare (= pulse) gives rise to 
a "pulse avalanche", and the whole light curve is produced by an 
overlap of such avalanches. Stimulated pulses are delayed from the 
stimulating pulse by a time At which is drawn from a Poisson dis- 
tribution of mean ari , where a is a delay parameter and n is the 
time constant of the stimulating pulse. This could be interpreted as 
a "reservoir" effect (Begelman & de Kool 1991). A long (i.e. en- 
ergetic) flare would drain the system of a larger amount of stored 
energy, which results in a longer time to store up enough energy 
for a subsequent flare. The presence of such reservoirs is also sup- 
ported by the observations (Negoro et al. 1995). The time constant 
of energy dissipation, r, for all pulses (spontaneous and stimulated) 
is assumed to be distributed according to a power law, p(r) oc r~ p 
between r m i n and r max . We assume that all pulses have the same 
shape and amplitude as shown in Figure 1 and differ only in time 
constant, r. In total we have six parameters: A, pi, a, r m i n , r max 
and p, but it is quite easy to navigate in the parameter space. In or- 
der to describe the role of the parameters, we discuss how each of 
them influence the PDS of a simulated light curve. 



5 RESULTS OF SIMULATIONS 

We assume that all flares produce the time-resolved spectra and 
energy-dependent pulses shown in Figures la,b. The variability 
and interaction between individual flares are described by a pulse 
avalanche model. Thus we are now able to simulate the light curves 
(Fig. 2a). The model parameters (A = 45; pi = 0.8; a = 7; p = 1; 
T m in = 1ms; and r max = 0.15 s) were selected to match a typ- 
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ical PDS of Cyg X-l, which can be described as a doubly-broken 
power-law: a flat, white noise, part f° below some frequency /flat; 
flicker noise, f 1 , between /fl a t and /t, r ; and "red" noise, / -/3 
with /3 ~ 1.5 - 2, above / br . Normally, / flat ~ 0.05 - 0.1 Hz and 
/ br ~ 2 Hz (Fig. 2b; Belloni & Hasinger 1990; Nowak et al. 1999). 

We identify /b r with the longest time-scale of an individ- 
ual flare r max « l/(27r/b r )- We note that flares of duration r, 
with the dissipation varying as (t/r)* exp(— t/r), produce PDS 
oc [(27r/) 2 + (l/r) 2 ]- ( « +1) decaying as / -2 « +1 ) at high frequen- 
cies and giving contributions to the total PDS in a very narrow fre- 
quency interval around / ~ 1/ (27rr) > / br , so that an actual value 
of £ affects none of the qualitative results. The absence of a break at 
/ ~ l/(27TT m j n ), corresponding to the minimum time scale, in the 
Cyg X-l data, constrains r m in ^> 2 ms. The flare time-scales can be 
related to the Keplerian time-scales at radii between the innermost 
stable orbit and ~ 30R g . The low frequency break /fl at ~ 0.1 Hz 
is related to the duration of the avalanche: ~ oiffj,/ (1 — /i) ~ Is, 
where f is the averaged over p(r) time-scale of the flare. The 
power-law index p of the time-scale distribution is related to the 
slope, (3 — 3 — p, of the PDS at / > /b r (e.g., Lochner, Swank & 
Szymkowiak 1991). The hardening of the PDS with photon energy 
(Nowak et al. 1999) is not reproduced by our model which predicts 
the same /3 for all energies. 

Parameters a and fi determine the duration of the avalanches, 
which, in turn, define the length and slope of the flicker noise part 
of the PDS, where most of the power emerges. The 1// part of the 
PDS is present for a very broad range of parameters (4 a ^ 10 
and 0.6 \i ^ 0.9). The larger a and fi, the longer the avalanches, 
the smaller /flat, and the larger the fractional root-mean-square 
(rms) amplitude of the variability (i.e. the normalization of the 
PDS). Increase in both a and /x, thus can account for the observed 
anti-correlation between /flat and the rms amplitude (Belloni & 
Hasinger 1990). The average number of spontaneous pulses, A, is 
a trivial parameter which defines the rms amplitude, rms oc 1/ \f\ 
(~ 30% in the case of Cyg X-l, Belloni & Hasinger 1990; van 
der Klis 1995). The rms amplitude decreases slightly with photon 
energy in agreement with the data (Nowak et al. 1999). 

Besides the PDS, we also compute the normalized cross- 
correlation functions (CCFs), time/phase lags, and coherence func- 
tions. The CCFs (Fig. 2c) are slightly asymmetric (see, e.g., Nolan 
et al. 1981) and peak at almost zero lags. They are well described by 
a stretched exponential law exp(— |i/fo|") with the index v ~ 1/2 
and to ~ 0.4 s within \t\ ^ 2 s. A similar stretched exponential 
behaviour is found in the autocorrelation function of gamma-ray 
bursts (Stern & Svensson 1996) and is known to describe some 
characteristics of turbulence (see, e.g., Ching 1991). 

Time lags, St(f) (Fig. 2d), decay approximately as 1// above 
/b r where individual pulses dominate the PDS: 

»</>"««•{ fit (2) 

where 5t max ~ r 'g ax \sx{E/Eo) is the maximum achievable time 
lag (here E and Eq are the corresponding photon energies). Since 
time lags are proportional to the time scale of the pulse, r ~ 
1 / (2-7r/), which gives contribution to the PDS at frequency /, they 
are smaller at higher frequencies. The maximum lag <5t m ax is a 
factor of two smaller than expected in a shot noise model with a 
fixed time-constant r max (see Fig. lb) due to averaging over dif- 
ferent time-scales. A wide distribution of time-scales of the indi- 
vidual flares ensures that there are non-negligible time lags and, 
at the same time, CCFs peaking at almost zero lags. The logarith- 
mic dependence of time lags on E is consistent with the observa- 



tions of Cyg X-l and other GBHs (e.g., Cui 1999; Nowak et al. 
1999). An amusing fact is that a similar logarithmic relation holds 
for Comptonization in a static cloud (Payne 1980) where higher 
energy photons, which spend more time in the medium, lag behind 
softer ones. The data for Cyg X-l (Nowak et al. 1999) show more 
complicated behaviour of the time-lags with Fourier frequency. Our 
model reproduces successfully 1/f dependence above ~ 1 Hz and 
a "shoulder" at ~ 0.5 — 1 Hz. A break in the time lag at ~ 0.5 Hz 
is not reproduced by the model. In § 6, we speculate on a possible 
origin of such a break. 

The coherence function (cf. Bendat & Piersol 1986; Vaughan 
& Nowak 1997) is close to unity (Fig. 2e), since the light curves 
at different energies are almost perfectly synchronized. Deviations 
from unity become visible at frequencies where individual pulses 
contribute to the PDS (/ <; /br ~ 1 Hz) and large E which is 
consistent with the data for Cyg X-l at / ^ 10 Hz (Nowak et al. 
1999), but not as strong as might be seen for / 10 Hz (although 
we note those data are especially prone to systematic uncertainties 
that were not fully considered by Nowak et al.). The larger E, the 
smaller the coherence function and the smaller Fourier frequency 
where loss of coherence starts. The same behaviour is expected for 
any energy dependent shot noise model. 

The parameter £ of the individual pulse profile affects the time 
delays between photons of different energies. This, in turn, affects 
the time lags and the coherence function in a way that the larger f , 
the larger the coherence and the smaller the time lags. Changing a, 
fl, p, or A in the pulse avalanche model hardly has any influence 
on the dependence of the time lags and coherence function on fre- 
quency. We note that changing r ma x affects the break in the PDS, 
/br ~ l/(27TT m ax), and correspondently changes the frequency 
where loss of coherence occurs. The maximum achievable time lag 
is also affected, <5t m ax oc r max . 

The total averaged luminosity of the system, L — 
£h.max~Ae 2 /[2(l — jtt)] ergs -1 , determines the average peak lu- 
minosity of individual flares, Lh,mnx- We find Lh,ma* ~ 1.6 ■ 
10 36 ergs -1 taking the Cyg X-l luminosity L = 4 ■ 10 37 ergs -1 
(Gierlinski et al. 1997). The peak compactness of the flare is then 
'ft, max ~ 15(R/ Rg)' 1 . This confirms a posteriori our choice of 
^h.max in § [| There are Ar/(1 — p,) ~ 7 flares active at any given 
time. The total energy dissipated during the flare Ea ~ L% )Tnax f 
constrains the magnetic field in the flaring region to be of the order 
B ~ ^6£fl/i? 3 ~ 10 s G. 



6 DISCUSSION 

The model we present can be developed in numerous ways. We as- 
sumed that the spectral evolution occurs due to the inflation of the 
magnetic loop configuration. This is, of course, not a unique possi- 
bility. A number of other scenarios can produce similar behaviour 
(see Poutanen & Fabian 1999). Another assumption that the spec- 
trum is formed instantaneously is no longer valid for short flares, 
r^a few R/c (~ 1 ms in our case, with a more stringent limit for 
higher observed photon energy), since delays due to Compton up- 
scattering and production of reprocessed soft photons then have to 
be incorporated. The delays in the production of soft (reprocessed) 
photons can make the Comptonized spectrum hard in the beginning 
of the flare and soft later producing negative lags at / <; 50 Hz. 

In our model, the longer flares appear to be very smooth, 
while, in reality, they can consist of a train of correlated, much 
shorter, flares. This effect can also significantly affect the PDS be- 
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Figure 2. (a) An example of a simulated light curve at 3 keV. (b) The power-density spectrum (PDS) at 3 keV (solid curve) and 27 keV (dashed). For a very 
broad range of parameters, the PDS can be represented as a doubly-broken power-law. The slope (3 at high frequencies is related to the slope of the time 
constant distribution, (3 = 3 — p (in our case /3 = 2). (c) Normalized cross-correlation functions for 9 keV (solid curve) and 81 keV (dashed curve) vs. 3 
keV. The 9 keV CCF peaks at a lag < 10~ 3 s and is almost symmetric; for higher energies, the peak of the CCFs shifts more towards positive lags and the 
CCFs become more asymmetric. Insert shows CCF behaviour near the peak, (d) Time lags, <5t(/) = 8<p(f)/(27rf) (thick curves), and phase lags, S(j>(f) (thin 
curves), for various energy bands: 9 keV (solid curve), 27 keV (dotted), and 81 keV (dashed) vs. 3 keV. Maximum time lags can be approximately described 
by 5i max ~ r °y x \n(E/Eo), where E and Eo are the corresponding energies and r max is the maximum flare time-scale. Since time lags are proportional 
to the time-scale of the flare (see Fig. lb), <5t(/) oc r ~ l/(27r/) (i.e. phase lags are almost constant) at frequencies / > / br related to individual flares, 
(e) Coherence functions, with the same definitions as (d). Note the different behaviour at frequencies below and above f bl and the decrease of the frequency 
where loss of coherence occurs when E increases. 



haviour at high frequencies (possibly even accounting for the hard- 
ening of the PDS with photon energy, Nowak et al. 1999). 

Another simplification is that we assumed that Compton re- 
flection is produced in the vicinity of the flare. However, it is pos- 
sible that some fraction of the hard radiation is reflected, e.g., from 
the outer disc, the wind, or the companion star. This would cause 
a suppression of the variability of this component at short time- 
scales (acting as a low-pass filter) and produce additional time lags 
independent of Fourier frequency. A break in the time lag - Fourier 
frequency relation should then appear for energies where Comp- 
ton reflection is not negligible (e.g., comparing E <; 7keV and 
E Si 2 koV channels) and it can be used to determine the distance 
of the reflector from the X-ray source. We propose that the break 
in the time lags at / ~ 0.5 Hz observed in Cyg X-l (Nowak et al. 
1999) can be a consequence of this effect. 

An important issue is the small amplitude of Compton reflec- 
tion seen in Cyg X-l and other GBHs in their hard state (Ebi- 
sawa et al. 1996; Gierlihski et al. 1997; Zdziarski et al. 1998), 
Q./(2tt) w 0.3. These observations have been interpreted as an 
absence of reflecting material in the direct vicinity of the X-ray 



source (cf. Dove et al. 1997; Poutanen, Krolik & Ryde 1997). In the 
present model, we explicitly assume that the cold material exists in 
the vicinity of the flare and produces soft photons by reprocess- 
ing hard X-ray radiation. Possible solutions include the following 
alternatives: (a) the cold disc is disrupted into a number of cold fil- 
aments that produce enough seed soft photons for Comptonization, 
but have a too small covering factor to produce visible Compton re- 
flection (Lightman 1974; Krolik 1998; Zdziarski et al. 1998; Pouta- 
nen 1998); (b) there is bulk motion directed away from the disc and 
more flare radiation is beamed toward the observer than the disc 
(Beloborodov 1999); or (c) the disc is highly ionized (Ross, Fabian 
& Young 1999). 

An interesting question is why the time lags in the soft state 
of Cyg X-l are smaller than those in the hard state (Cui et al. 
1997). We can suggest at least two reasons. First, since the soft 
state spectrum is softer, the amplitude of the spectral evolution is 
much smaller. Second, if in the soft state X/7-rays are produced 
only within ~ 10R g from the black hole, the dynamical time-scales 
are on average smaller, causing smaller lags. 
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7 SUMMARY 

We have presented a model for the short time-scale spectral vari- 
ability of accreting black holes, and showed using Cygnus X-l that 
it is able to describe time-averaged spectra as well as general be- 
haviour of the temporal characteristics such as power-density spec- 
tra, time/phase lags, and coherence function. We based our spectral 
model on the fact that magnetic loops have a tendency to inflate and 
detach from the accretion disc. The corresponding decrease of the 
feedback from the cold disc causes the soft-to-hard evolution of the 
flare spectrum which is the reason for time lags between hard and 
soft photons. We assumed that all properties of the flares scale with 
their time-scale, r. This produces the time delays between hard and 
soft photons in an individual flare also scaling with r. In their turn, 
the time lags from an ensemble of flares of different duration show 
the 1// dependence on the Fourier frequency. 

The model for the temporal variability is based on a pulse 
avalanche model where flares can, with some probability, trigger 
other flares. The duration of the avalanches, which are the result of 
this interaction, determines the Fourier frequencies where most of 
the power emerges. 

We identified the break in the PDS at f^ t with the maximum 
time-scale of the individual flare r max « l/(27r/b r ) and pointed 
out that statistical properties are very different at / < /br and 
/ > /br- Deviation from unity of the coherence function and the 
break in the time/phase lag spectrum occur also close to /br- The 
PDS, time/phase lags, and coherence at / > /b r are related to the 
properties of the individual flares, while the behaviour at / < /b t 
is a property of a stochastic process. Our model predicts that vari- 
ations of /br in the PDS should be accompanied by corresponding 
changes in the time/phase lag and coherence spectra. 

The success of the model in reproducing most of the spectral 
and temporal properties of Cygnus X- 1 supports the point of view 
in which X/7-rays are produced in localized hot regions (magnetic 
flares) that interact with each other. We expect that this model may 
also be useful when analyzing data from other accreting objects 
such as active galactic nuclei and weakly magnetized neutron stars. 
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